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PREFACE 


Professor  Norman  D.Ham  is  presenting  here  a  revised  and  up-dated  version  of  the  article  he  wrote  in  1967, 
under  the  same  title,  for  the  Aeroelasticity  Manual,  and  which  was  included  in  Chapter  10  of  Volume  III  the 
following  year. 

Since  that  date,  many  advances  and  developments  have  been  made  regarding  the  vibration  theory  of  the 
rotating  parts  of  helicopters,  and  the  understanding  of  the  instabilities  created  in  flight  by  these  vibrations.  For 
the  Manual  to  continue  to  fulfil  its  information  mission,  it  had  become  necessary  to  bring  it  up-to-date.  Professor 
Ham,  who  had  been  the  first  one  to  draw  attention  to  developments  in  the  field  covered  by  his  article,  accepted 
this  task  himself;  no  one  could  have  brought  it  to  a  more  successful  issue. 

New  questions  i«a'v  been  raised  and  solved.  The  analytical  methods  used  to  convert  in-flight  vibrations  into 
equations  and  to  detect  aeroelastic  instabilities  are  presented  for  hinged  and  locked  blades,  for  take-off  or  cruise 
flights.  The  reasons  for  such  instabilities  are  considered  in  detail,  and  efficient  means  prescribed  to  avoid  them. 
Besides  conventional  type  flutter,  a  large  section  is  devoted  to  stall  flutter.  The  text  is  illustrated  and  completed 
by  many  diagrams  giving  the  results  of  calculations  carried  out  in  the  United  States. 

Considered  from  a  general  standpoint,  this  article  which  is  intended  to  replace  that  of  1968,  provides  an 
excellent  survey  of  the  data  available  in  1972  on  the  vibratory  stability  in  flight  of  hel  copters,  and  should  therefore 
prove  extremely  useful  to  helicopter  engineers  and  manufacturers. 


R.  MAZET 
General  Editor  of  the 
Manual  on  Aeroelasticity 


PREFACE 


Le  Professeur  Norman  D.Ham  presente  l:i  une  edition  revisee  et  compl6tee  de  Particle  redige  par  lui  sous  le 
me  me  ti  re  en  1967  pour  le  Manuel  d’Aeroelasticite  et  qui  a  pris  place  au  Chapitre  10  du  Volume  III  I’arinee 
suivante. 

Depuis  cette  date,  de  nombreux  progres  et  developpetnents  avaient  ete  apportes  a  la  theorie  des  vibrations  des 
organes  toumants  des  helicopteres  et  a  la  connaissance  des  instabilities  auxquclles  ces  vibrations  donnent  naissance 
en  vol.  Pour  que  le  Manuel  continue  de  remplir  la  mission  d'information  qui  lui  a  ete  assignee,  il  devenait  necessaire 
de  proceder  a  une  remise  a  jour.  Le  Professeur  Ham.  qui  avail  le  premier  signale  1’evolution  des  connaissanccs  sur 
la  matiere  de  son  article,  a  bien  voulu  assumer  lui-meme  ce  travail:  nul  ne  pouvait  mieux  que  lui  le  mener  a  bien. 

Des  questions  nouvelles  sont  posees  et  residues.  I.es  me'.hodes  analytiques  permettant  de  mettre  en  equations 
les  vibrations  en  vol  et  de  deader  les  instabilities  aeroelastiquvs  sont  presentees  pour  des  pales  articulees  ou 
encastrecs.  pour  le  vol  au  decollage  ou  en  croisiere.  Les  raisons  de  ces  instabilities  sont  examinees  en  detail  et  des 
moyens  efficaces  sont  preconises  pour  les  eviter.  A  cole  des  flottements  de  types  classiqucs.  une  large  place  est 
fait  au  ffottement  de  deerochage  (stall  flutter).  De  nombreux  graphiques  traduisant  en  courbes  les  resultats  de  calculs 
effectues  aux  USA  illustrent  et  completent  le  texte 

D’une  fa  yon  generate,  eet  article,  qui  est  destine  a  se  substituer  a  celui  de  1968.  est  une  exeellente  synthese 
des  connaissances  disponibles  en  1972  sur  la  stahilite  vibratoire  en  vol  des  helicopteres  et  doit  rendre.  a  ce  litre, 
d’utiles  services  aux  const ructeurs  et  aux  ingilnieurs  splcialistes  de  ces  appareils. 


R  MA7I.T 

l.diteur  General  du  Manuel 
d'Aeroelasticitc1 
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HELICOPTER  BLADE  FLUTTER 

by 

Herman  D.  Ham 
Professor 
MIT 

Cambridge,  Mass.  OT139 
U  .  u .  A . 


BIMHaRY 


Methods  of  analysis  of  nellcopter  blade  flutter  for  both  hinged  and  hingeless 
blades  are  presented.  The  major  types  considered  are  bending-torsion  flutter,  flap-lag 
flutter,  and  stall  flutter.  Both  tiover  and  fnrwaru  flight  are  considered.  Means  of 
avoiding  flutter  are  described. 


SYMBOLS 
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g(t ) 
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m(r) 
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r  or  s 
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M(r) 
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n(  r) 
nk(r) 
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blade  section  lift  curve  slope,  per  radian 

blade  section  seni-cnord 

blade  section  chord 

displacement  of  blade  first  elastic  bending  mode 

displacement  of  blade  Kth  bending  mode 

blade  section  reduced  frequency,  - — - 

i2R(x  +  u  slrul) 

blade  spanwise  running  mass,  slugs  per  foot 
blade  nondlmenslonal  spanwise  station 

cnorawise  displacement  of  section  aerodynamic  center 
from  elastic  axis,  positive  forward 

cnorawise  displacement  of  section  center  of  gravity 
from  elastic  axis,  positive  forward 

spanwise  distance  along  blade  from  axis  of  rotation 

vertical  distance  from  rotor  nub  plane 

lift  deficiency  function  Including  blade  wake  effects 

rotor  tnrust  coefficient 

blade  bending  modulus  of  elasticity 

real  part  of  lift  deficiency  function 

Imaginary  part  of  lift  deficiency  function 

blade  section  second  moment  of  area,  flatwise  bending, 
or  ratio  Ig/I^ 

blade  moment  of  Inertia  about  flapping  ninre 

blade  section  moment  of  Inertia  about  its  center  of  gravity 

nondlmenslonal  blade  product  of  inertia  auout  flapping  anu 
elastic  axes 

nondlmenslonal  generalized  mass  of  kt.n  benuing  node 
blade  section  moment  of  Inertia  about  its  elastic  axis 
blade  bending  moment  at  r 

blade  section  aerodynamic  moment  about  elastic  axis 


rotor  radius 
rotor  thru?’ 

blade  section  angle  of  attack 

blade  section  mean  angle  of  attack 

displacement  of  blade  rigid  flapping  mode 

*i 

blade  mass  constant,  or  Lock  number,  - 

*1 

mode  shape  of  blade  first  elastic  bending  node 
mode  shape  of  blade  kth  bending  mode 

rotor  advance  ratio  V/i2R 
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v  nondimenslonal  root  of  bending-torsion  characteristic  equation 

nondimenslonal  rotating  natural  frequency,  kth  bending  mode 
u  imaginary  part  of  v 

rotating  natural  frequency,  kth  bending  mode 
nonrotating  natural  frequency  of  blade  torsional  motion 

p  air  density 

o  rotor  solidity 

0  blade  pitch  angle 

0  damping  coefficient  of  blade  torsional  motion,  divided  by 

tne  critical  damping  coefficient 

Si  rotor  rotational  speed,  radians  per  second 

iji  rotor  blade  azimuth  angle,  zero  when  reference  blade  is  downstream 

©  (x)  fundamental  torsion  mode  shape 


1.  INTRODUCTION 

Tne  flutter  of  helicopter  blades  can  be  classified  under  the  following  categories: 

(I)  Flutter  of  Hinged  Blades 

(J)  Flutter  of  Hingeless  Blades 

(3)  Stall  Flutter 

The  analysis  of  case  (1)  differs  from  tne  classical  analysis  of  the  high  aspect 
ratio  wing  due  to  several  important  factors.  The  high  centrifugal  force  field  experienced 
by  tne  rotating  blade  provides  fundamental  bending  and  torsional  stiffness  effects  ?.nd 
also  important  coupling  effects  under  certain  circumstances.  Also,  the  presence  of  the 
retaining  wake  beneatn  the  rotor  at  low  flight  speeds  leads  to  unsteady  aerodynamic  ef¬ 
fects  tnat  are  substantially  different  from  those  characteristic  of  fixed  wings.  Finally, 
tne  variable  velocity  field  encountered  by  the  blade  in  forward  flight  generates  Important 
time  variations  in  the  velocity-  and  amplitude-dependent  restoring  forces  acting  on  the 
blade . 


The  flutter  In  case  (2)  is  due  to  the  coupling  between  blade  flapwise  and  lagwise 
motion  whic.i  may  lead  to  instability  at  large  blade-pitch  settings.  This  Instability  can 
occur  at  any  flight  speed  provided  the  damping  of  the  chordwise  motion  is  below  a  certain 
value . 


The  instability  associated  with  case  (3)  Is  due  to  the  adverse  time  phasing  of 
the  aerodynamic  torsional  moment  resulting  from  the  loss  of  blade  bound  vorticity  during 
torsional  motion  at  nign  angles  of  attack.  The  complex  nature  of  the  phenomenon  precludes 
an  analytical  representation  of  the  unsteady  airloads  at  the  present  time,  but  the  pre- 
uiction  of  regions  of  instability  is  made  possible  by  the  application  of  experimental  two- 
dimensional  data. 

A  useful  survey  of  tne  various  instabilities  of  rotors  is  presented  in  Refer¬ 
ence  ( 1)  . 


2.  FbuTTLR  OF  HINGED  BLADEJ 

a.  Equations  of  Hotlon 

Consider  bending  out  of  the  plane  of  rotation  of  the  flexible,  untwisted,  un¬ 
tapered  rotating  blade  snown  In  Figure  (1).  The  bending  moment  at  r  due  to  forces  at  s 
is 


M  ( r ) 


d  T  ( s ) 
ds 


-  z  ( S )  m  ( s ) 


(s-r) ds 


m  (s) 


z 


ds 


differentiating  twice  with  respect  to  r  yields 


d2M!r) 
d  r 2 


d  T(r ) 


m  ( r  )  7.  ( r )  + 


d  2  Z  ( r ) 


m (s) ds  - 


dz(  r ) 


dr 


dr 


ril2  m  ( r) 


r 

i 


M (r)  =  E I (r) 


d2z 

dr2 


dr* 


El  (r) 


d2z 

dr2 


-ill  f 

dr2  J, 


m  (s)  sfl2  ds  +  rfl2  m(r)  —  +  m  ( r)  z  = 


dr 


d  r 


(1) 


the  differential  equation  for  a  beam  bending  in  a  centrifugal  force  field. 
Assume  a  series  solution  in  terms  of  normal  bending  modes. 


co 

2  -  Z  T)  (r)  g  (t) 

k  =  l 


(2) 


For  free  vibration  of  the  rotating  beam,  neglecting  aerodynamic  damping, 
dT(r)/dr  «  0.  Substituting  Eq.  (2)  into  Eq.  (1) 
j2_  j2_  a  R 


r° 


quency 


Assuming  simple  harmonic  motion  at  the  kth  rotating  undamped  natural  bending  fre- 

_  ii/kftt 

9k  8  9k  e 


v _  at 


-"k2ft29k 


[  1 


vl  fl2m,?k  9k 


Therefore,  the  bending  equation  of  motion  becomes 


oo 


v  2  n2  S  .  dT 

>  mu  D  g  t  )  rrrn,  g,  =  -7— 
k  'k  *k  .  ,  'll  k  dr 
k:|  k-l 


As  an  approximation,  assume  that  all  blade  torsional  flexibility  is  concentrated 
at  the  root  (often  true  due  to  control  system  flexibility:  the  general  case  of  distributed 
torsion  is  discussed  in  Reference  (2)).  Then,  for  the  geometry  of  Figure  (2),  the  coupled 
blade  bending-torsion  equation  of  motion  is 


CO 

z 

k=l 


CO  c  c 

Z  m  17  a  +  Z  mi/  A  17  a  -mx.0-mx.fl  8  -  -3 — 
7,  'k  Vk  k71  k  'kvk  I  l  dr 


2  r.2 


dT 


(3) 


Multiply  all  terms  of  Eq,  (3)  by  and  Integrate  from  0  to  R: 


Mk  9'k  +  Mkyk 


fl2gk-  (9  +  fl2ff  )  J  rn-»7k  *  j  d  r  -  J  dr 


(4) 


where 


M  «  /^  2 

nk  'o  dr 

/; mn  j  nv< dr  *  ° 


k  -  1,  2,  .. 
J  ¥  k 


due  to  the  orthogonality  of  the  n^'s. 

In  practice  an  adequate  representation  of  the  blade  bending  motion  is  often  ob¬ 
tained  in  terms  of  the  first  two  modes  only.  On  this  basis  the  bending  motion  of  an 
articulated  blade  is  described  by 


z  -  r0  +  ng 


(5) 


where  r  ■  rigid  flapping  mode 
n  "  first  elastic  mode 


4 


Then  if  I  s  ;  I2  ■  M^/I  ;  Ix  ”  (l/I^)/^  mn^Xj  dr»  Ec's*  (**)  become,  after  division  by 
1^  Si2,  and  since  *  1, 


0_ 

ft2 


+  0 


(6) 


IT 


hvh 


- 1, 


'2  n‘ 


- 1 


*2e  = 


rR  di 

J  v^~ 


dr 


(7) 


How  consider  torsional  motion  of  the  blade  shown  in  Figure  (2).  Taking  moments 
about  the  elastic  axis, 


-I 


ft  R  R 

m  (z-Xjfl)  x  j  dr  +  f  IQSdr+  f  lo&20dr  +  I  w  2  0 
0  *0  Jo 


fR  n2  d z (o)  fR  n2  dz(o)  .  l  /*R  dMft, 

-J  rmn  ~dr~  XI  dr  ~J  mxrn  (2"r-d7--xi  6]dr'-n?J  ~d 


I  0 


n  p  0 

for  small  z  and  Q,  anu  no  lag  hinge.  The  term  [I  +  f  xT  dm]  Si  6  is  the  well-known 
"propeller  mo  went". 

Substituting  Eq.  ( b  >  and  dividing  by  i^Sl2, 


(8) 


The  aerodynamic  terms  in  Eqs.  (6),  (7),  and  (8)  are  now  considered  in  detail. 

From  Reference  (j),  P.  272,  for  unU  span  of  a  thin  airfoil  oscillating  in  incom 
pressible  flow,  and  neglecting  small  unsteady  effects  due  t-o  tne  unsteady  component  of  U_ 
Reference  (4), 


dj 

dr 


-  -  y  pac2  [z  +  UjS  -  ( xft-  0.25  c  )  8  ] 

-i  pocUTC'(k)  [up  +U j6  +  (0.5c-xa)  e] 


dM*  i  ?  r  •  c2  •  2  i 

— —  -  —  pa  cZ  [(  xA-  0.2  5  c  )  z  -  UT  9  (  0.5  c-  xA)-  yg-  0  -  (xft-  0.25c  )  0  j 


p  ac  xft  C*  (k)  [up+UT0  +  (0.5c  -  xft  )  0  ] 


where  the  blade  incident  free-stream  velocity,  including  rotor  rotational  velocity,  is 

U,^,  *  Sir  +  pSiR  siniii 

and  the  unsteady  flow  perpendicular  to  the  blade  is 

•  d  z 

Up  *  z  +  pSiR  cosi^ 


The  quantity  C'(k)  is  comparable  to  the  classical  Theodorsen  function  C(k),  but  has  im¬ 
portant  differences  due  to  the  unique  characteristics  of  the  rotor  wake.  Typical  values 


5 


of  C'(k)  are  shown  in  Figure  (3),  taken  from  Reference  (5). 

In  the  present  analysis,  the  virtual  mass  terms  proportional  to  S !  and  0  are 
neglected  for  simplicity.  They  can  be  included  if  desired  by  appropriate  adjustment 
(usually  of  the  order  of  a  few  ner  cent)  of  the  values  of  the  blade  inertial  constants. 

The  aerodynamic  terms  can  be  expressed  as  follows: 


1  fR  dT  .  P  Q  on 

i  a2  J  r  dr  dr  '  "  m/j/n  si  m2  si  "  me/n  ”  me  9  "  m0  ^  '  m|*2  9 


where 


m . 


t 5/Xl  8 


m 


m 


m 


m 


m 


=  “  C'(k) 


f1  f  H-  S'"'/'] 


*  Y  C'(k)  [  J  %■  dx  +  fj  sin  \p  J  *  d x  J 

1  £[,  +  7*sH[i  +  4  (05f-ir)  ^>] 

s^c'(k)  [(  ’  +  ^.2)  +  -|-  /isin<//  -  /x2cos  2<|r  ] 

!  y  C'(k)  [y  /i  cos  if/  +  sin  2  ip  j 

*  f  C'(k)  [MCOS^/lx2^r(^-)dx  +  ^M2sin2^/,x^r 


The  Q/Q  term  is  normally  neglected  for  c  <<  R.  The  coefficient  C' (k)  is  a  mean  value 
based  on  condition?  at  the  blade  three-quarter  radius.  Also 


/•  R  "  •  * 

_L_  |  v  dldrs.m2  •f-^20-m^/3-mgg 


where 


Vfl 


j/«  s  f  c,,kl  ♦Msr+jf'dfd-] 

'  -f  [m *i" *•][-§■’- «  (0.5|-^l  + 

=  Y  c'(k)  [  J*  (Aj  x2dxr^2sin2^J*  -j|dxj+yC'(k)/isin</rJ*  x  dx 
1  Y  C'(k)  ^  nZ  sin  2  ^  J*  -2-dx  +  ^  cos  ^  *  -R-dxj 

-  £■  C'(k)  [mcos^  ,^A^dx  +  ±^sinZf^  (-i)A(i)dx] 


Again,  the  0/n  term  is  neglected  for  c  <<  R.  Finally, 
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Note  that  in  forward  flight,  when  the  flow  direction  reverses  periodically  with 
respect  to  portions  of  the  blade,  x,/R  ■  0.5  c/R,  and  the  blade  damping  in  pitch  Mi 
becomes  zero  for  such  portions . 

b.  Method  of  Analysis 
(i)  In  Hover 


For  discussion  purposes,  torsional  motion  0  and  rigid  flapping  motion  6  of  the 
rotating  blade  will  be  considered.  Then  the  equations  of  motion  become 


*  m/j/a  Q  +  (1  0  -  ♦  <m,  _I«  >  ®  ‘  0 

[V'*‘  (f1)2]®10 


In  hover,  these  equations  have  constant  coefficients  and  can  be  solved  by  con¬ 
ventional  techniques.  Assuming  simple  harmonic  motion,  where 


B  -  ffevnt 
0  -  0evIJt 

the  characteristic  equation  of  the  motion  is 

Av^  +  Bv3  +  Cv2  +  Dv  +  E  ■  0 

For  the  cace  of  quasi-static  airloading,  C'  (k)  *  1,  all  coefficients  of  the  characteristic 
equation  are  real,  and  the  stability  of  disturbed  motion  can  be  evaluated  using  Routh's 
criteria. 


When  rotor  wake  effects  are  Included,  C' Ik)  is  complex  and  a  function  of  v,  and 
a  trial  anc  error  solution  is  necessary.  For  the  case  of  neutral  stability,  v  *  iw, 
where  u  is  tne  flutter  frequency,  and  the  characteristic  equation  can  be  separated  into 
its  real  and  imaginary  components, 

A'w1*  +  B'w3  +  Cw2  +  D'w  +  E'  •  0 
A'  '(i)1*  +  B''u3  +  C’u2  +  D' 'w  +  E"  -  0 

where  the  coefficients  are  functions  of  u.  The  pair  of  equations  can  then  be  solved  simul¬ 
taneously  by  assuming  values  of  u  and  solving  each  of  the  two  equations  for  the  parameter 
2 

(«0/Q)  .  Values  of  u  which  yield  identical  values  of  this  parameter  for  the  two  equations 
are  the  flutter  frequencies  Up,  and  the  torsional  stiffness  or  rotational  speed  at  flutter 
is  defined  by  the  corresponding  value  of  uQ/n.  Conditions  for  torsional  divergence  are 
not  dependent  on  wake  effects  since  for  divergence  u  -  0. 

The  flutter  and  divergence  boundaries  for  a  hovering  rotor  were  computed  in 
Reference  (6)  for  a  rotor  blade  having  the  following  properties: 

C’TTc)  1 

Y  -  12 
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The  results  are  presented  In  Figure  (A). 
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The  effect  of  the  wake  Is  shown  in  Figure  (5)  using  values  of  C'(k)  from  Refer¬ 
ence  (7)  to  recompute  the  hovering  boundary  of  Figure  (4).  It  is  seen  that  neglecting 
the  effect  of  the  shed  wake  Is  adequate  for  design  purposes  in  this  instance. 


The  effect  of  blade  flexibility  on  the  flutter  boundaries  i3  also  shown  in 
Figure  (5)  by  a  single  point  computed  for  a  blade  bending  frequency  v,  ■  2.5-  For  con¬ 
ventional  blades,  neglect  of  bending  flexibility  results  in  negligible,  though  unconserva- 
tlve,  error. 


Further  discussion  and  comparisons  with  test  results  are  presented  in  References 
(8)  and  (9). 


Usually  flutter  of  conventional  articulated  helicopter  blades  is  avoided  by  mass 
balancing  them  about  the  quarter-chord  point,  i  e. ,  by  making  the  blade  chordwiLe  center- 
of-gravity  position  coincident  with  the  blade  section  aerodynamic  center.  However,  in 
certain  Instances  appreciable  blade  bending  coupled  with  a  forward  elastic  axis  location 
can  lead  to  flutter  even  with  mass  balanced  blades,  References  (8)  and  (10). 

(li)  In  Forward  Flight 

The  equations  of  Section  2a  can  be  solved  numerically  for  the  actual  blade  motion, 
including  the  forcing  aerodynamic  terms  due  to  blade  pitch  angle,  variable  downwash .steady 
state  flapping  and  bending,  advancing  blade  Mach  number,  stall,  and  reverse  flow. 

A  recent  Investigation,  Reference  11,  presents  step-by-step  computations  of  rotor 
blade  motion,  Including  the  effect  of  blade  dynamic  stall,  during  flight  at  moderate  to 
high  advance  ratios.  The  three  degrees  of  freedom  considered  in  the  analysis  are  rigid 
blade  pitching  about  the  feathering  axis,  rigid  blade  flapping  about  the  zero  offset 
flapping  hinge,  and  blade  first  mode  flatwise  elastic  bending.  Variable  Inflow  and  re¬ 
verse  flow  are  included  in  the  computation  of  blade  forces  and  moments.  The  calculation 
of  Inflow,  aerodynamic  loading,  and  blade  motion  is  performed  iteratively,  using  estimated 
steady  state  initial  values. 

Typical  results  are  shown  in  Figures  (6)  and  (7).  The  nature  or  these  results 
suggests  that  classical  methods  of  blade  flutter  and  divergence  analysis,  carried  over 
from  fixed  wing  practice,  should  be  replaced  by  numerical  solution  of  the  blade  equations 
of  motion  including  not  only  the  conventional  harmonic  excitation  forces,  but  also  the 
nonlinear  effects  of  reverse  flow,  variable  Inflow,  dynamic  stall,  and  compressibility. 

Classical  flutter  is  a  highly  periodic  self-excited  motion;  rotor  blade  motion 
at  high  advance  ratio  is  largely  due  to  external  excitation.  The  actual  level  of  blade 
motion  is  the  important  criterion,  as  seen  in  the  present  results,  which  demonstrate  not 
a  stability  boundary,  but  a  blade  stress  level  boundar"  due  to  the  excited  motion.  Even 
the  "classical  flutter"  case  exhibits  the  effects  of  first  harmonic  excitation,  while  the 
"torsional  divergence"  case  is  the  blade  response  to  the  lift  acting  at  the  blade  three- 
quarter  chord,  and  varying  with  the  dynamic  pressure  over  the  reverse  flow  region. 

At  very  high  advance  ratios,  uncoupled  flapping  Instability  can  occur  due  to 
negative  aerodynamic  damping  and  spring  forces  resulting  from  the  periodic  nature  of  the 
flow  (Reference  12). 


c.  Finite  Bending  Effects 

The  preceding  discussion  has  neglected  the  effect  of  steady  state  coning,  collec¬ 
tive  pitch,  and  elastic  bending  displacement  of  the  blade.  For  a  conventional  articulated 
rotor,  these  effects  are  small.  However,  for  certain  hub  geometries,  the  steady  3tate 
displacements  of  significant  modes  are  of  fundamental  importance  and  cannot  be  neglected. 

A  detailed  analysis  of  such  effects  is  presented  in  Reference  (2).  For  example,  in  the 
absence  of  a  lag  hinge,  steady  state  flatwise  bending  displacement  of  the  blade  intro¬ 
duces  destabilizing  moments  due  to  centrifugal  force  which  in  extreme  cases  lead  to  in¬ 
stability  when  the  blade  chordwise  center  of  gravity  coincides  with,  or  is  even  ahead  of 
the  blade  aerodynamic  center.  These  destabilizing  moments  are  largely  attenuated  by  the 
incorporation  of  a  lag  hinge. 


3.  FLUTTER  OF  HINGELESS  BLADES 
a.  Equations  of  Motion 

Classical  blade  flutter  occurs  due  to  the  coupling  of  the  torsional  and  flapping 
degrees  of  freedom.  Under  certain  conditions,  the  two  bending  degrees  of  freedom  of  the 
blade  (in  the  plane  of  rotation  and  out  of  the  plane  of  rotation)  couple  together  to  pro¬ 
duce  another  type  of  instability  called  flap-lag  flutter. 

In  this  section,  the  flap-lag-type  Instability  of  torsionally-rigid  hingeless 
blades  in  the  linear  range  of  blade  motion  is  treated  in  hover.  This  problem  wp.i  first 
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treated  by  Young  with  a  restrictive  analytical  approach  (Reference  13) •  Modal  equations 
of  motion  were  obtained,  but  the  numerical  results  were  evaluated  for  a  blade  represented 
by  a  centrally-hinged,  spring-restrained,  equivalent  model.  Young  concluded  that  the 
triggering  mechanism  of  the  flap-lag-type  Instability  Is  the  lag  degree  of  freedom. 


Hohenemser  treated  the  same  problem,  using  a  somewhat  unconventional  numerical 
integration  scheme  (Reference  l**).  Due  to  the  various  approximations  made  In  Reference  H, 
the  results  presented  there  are  of  n  qualitative  nature. 

Consider  the  hingeless  bl».dt.  shown  in  the  hub  plane  axis  system  of  Figure  (8). 
Following  the  analysis  of  Reference  15,  the  displacements  v  and  w  are  expressed  as 


v/R 


-  y1(x)h1(t) 


w/R  ■  n1(x)g1(t) 


(1) 

(2) 


where  n^x)  Is  the  first  normal  mode  shape  of  blade  lapwise  bending,  ana  g^t;  is  the 
modal  displacement.  Similarly,  Is  the  first  normal  mode  shape  for  blade  lagwlse 

bending,  and  h^(t)  Is  the  modal  displacement. 

The  nonlinear  equation  of  motion  for  flapwise  bending  Is  snown  In  Reference  15 


to  be 


=  *  i  [F,  V  %  3,-|2F,e8-F„  Vh'  + 

i,*,] 
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to  be 


X0  =  rotor  inflow  ratio  referred  fo  hub  plone 
(*)  =  d/dVr  (  ) 1  =  d  /dx 

The  nonlinear  equation  of  motion  for  lagwlse  bending  Is  shown  In  Reference  15 


ML| ( TT(  +  w*  h,)*  (2S,  -  2  My(  )  h |  h  !  -  2M^  g,  g, 

+  \  [(L,e-L2  X0)  X04  -p-  L4 
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~  LI8  9I  ~  LI9®  9I  hl  ] 


where 


w, 


d  x 


first  mode  rotating  logging  frequency 


9 


§,  * 

V 

>.  Jl  V 

R3  /*' 

[/'-  X.  **,' 

[jN>v 

di 

1 

M  -- 

r, 

r3  r' 

T,  l  " 

di 

cd  * 
do 

blode  section 

profile  drag  coefficient 

L.  • 

r !  *  ai 

i 

Y  2  d  * 

L2 

r'r  di 

1 

^■r 

■  y,  d* 

L4  * 

/»'  _  2 
i  ’  r> 

l,b  'l 

V  r. d' 

L7  ‘ 

•'o 

L|9  '•( 

i?,  y 2  d* 

L8  * 

f'v.r,  ai 

Equations  (3)  and  (9)  are  a  system  of  two  second  order  nonlinear  equations  de¬ 
scribing  the  motion  of  the  syatem.  These  equations  are  coupled  through  the  following 
effects: 

(a)  Coriolis  forces 

(b)  Shortening  effects 

(c)  Aerodynemlc  forces 


The  nonlinear  system  given  by  Eqs.  (3)  and  (9)  are  linearised  about  the  static 
equilibrium  condition  (l.e.,  with  all  time  derivatives  set  equal  to  zero),  denoted  by 
gj  ,  hj  ;  from  the  eq.stions  It  Is  clear  that 
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Substitution  of  Eqs.  (7)  and  (8)  Into  (3)  and  (9),  and  use  of  Eqs.  (5)  and  (6)  yields, 
aft*r  neglecting  terms  of  type  Ag^  -Ah^,  etc.,  and  writing  Agj  ,  as  g^  for  convenience, 


M  (  g,  9,)-2  P,g 


ml,  hi  >s  -  2  Mt>,  9i0  St* 
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Note  that  the  term  2(?1  -  Hy^h^h.  ■  0  has  been  omitted  In  Eq.  (10)  because  It  Is  usually 


zero.  For  convenience,  the  following  terms  will  be  defined: 
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With  Eqs .  (11)  through  (14),  Eqs.  (9)  and  (10)  can  be  rewritten  in  the  following  conveni¬ 
ent  manner 
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Prom  Eq.  (15),  It  Is  clear  that  X,Y  represent  the  coupling  terms,  and  from  Eqs.  (12)  and 
(14)  it  Is  clear  that  the  coupling  Is  partly  due  to  Coriolis  effects  and  partly  due  to 
aerodynamic  effects. 


The  quantities  D^,  represent  the  damping  in  the  system,  the  damping  lnflap  Is 

a  relatively  large  number,  while  the  damping  tn  lag  Is  a  small  quantity  of  order  C.  /a; 
the  lag  degree  of  freedom  Is  the  potentially  unstable  one.  ao 

b.  Method  of  Analysis 

The  flutter  or  the  critical  condition  of  the  linearized  system  (15)  Is  character¬ 
ized  by  the  existence  of  a  small  amplitude  oscillation  for  Eqs.  (15). 

Assume  the  solution  In  the  form 


Ale 


P9 


*2e 


P* 


Substitution  of  these  Into  (15)  yields  the  following  characteristic  equation 


( P2  ♦  0,  p  ♦  )(p2  ♦  D2  p  -  p2  XY  *  0  (16) 

For  a  small  vslue  of  6,  the  root  of  Eq .  (16)  has  Real(p)  <  0  end  the  solution  Is  stable. 
Por  9  •  9c,  the  system  Is  neutrally  stable.  Por  0  >  9c,  at  lesst  one  of  the  roots  of 

Eq.  (16)  has  Real(p)  >  0  and  the  system  Is  unstable. 

At  9  ■  9c,  the-e  are  two  solutions  to  Eq.  (16)  such  that  p  Is  Imaginary: 


yields 


p  ■  iiw 

Substituting  Into  (16)  and  setting  to  zero  the  real 
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and  Imaginary  part  of  (16) 

(17) 

y  »0  (lfl) 


where  w  Is  the  flutter  frequency.  It  Is  Interesting  to  note  that  3-,/D,  <<  1.  Therefore, 

C  t  1 

2  v  -2 
“c  '  “LI 

Stability  boundaries  resulting  from  the  solution  of  Eqs.  (17)  and  (18)  are  shown 
In  Pigure  9  for  no  structural  damping  and  the  mode  shapes 

^  (1  -  4?  -  (1  -  a)*] 


The  areas  Inside  the  elliptical  boundaries  are  cos&lnatlons  of  rotating  flap  and 
lag  frequencies  for  which  the  system  Is  unstable.  These  areas  are  reduced  for  Increasing 
structural  damping,  and  Increased  for  Increasing  y. 

c.  Torsional  Effects 

The  addition  of  the  torsional  degree  of  freedom  Is  stabilizing  for  the  lower 
branch  or  the  flap-lag  stability  boundary  shown  In  Pigure  9,  and  strongly  destabilizing 
for  the  upper  branch,  as  shown  In  Pigure  10. 

Note:  Important  effects  are  associated  with  the  nonllnearltles  contained  in  Eqs.  (3)  and 
(4).  Por  further  details,  see  Reference  (15). 
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H .  STALL  FLUTTER 

a.  Aerodynamic  Loading  During  Dynamic  Stall 

Studies  described  In  References  (20)  and  (21)  have  shown  that  the  negative  damping 
In  pitch  associated  with  airfoils  oscillating  at  high  mean  angles  of  attack  can  lead  to 
torsional  Instability  of  helicopter  rotor  blades  under  certain  conditions.  The  mechanism 
of  the  instability  was  shown  to  consist  of  the  adverse  time  phasing  of  the  aerodynamic 
pitching  moment  associated  with  the  loss  of  blade  bound  vortlclty  as  the  dynamic  stall 
occurs . 


Subsequent  tests,  Reference  (22),  Indicated  that  the  same  mechanism  Is  found  for 
the  case  of  airfoil  linear  angle-of-attack  change  through  high  angles  of  attack.  The 
nature  of  these  results  Indicated  several  Important  conclusions  not  only  with  respect  to 
the  analysis  of  rotor  blade  stall  flutter  Instability,  but  also  with  respect  to  the  gen¬ 
eral  nature  of  the  aerodynamic  loading  of  an  airfoil  experiencing  transient  angle-of- 
attack  changes  of  large  magnitude. 

A  typical  time  history  of  the  pressure  variation  acting  at  one  spanwlse  station 
of  a  model  helicopter  rotor  blade  experiencing  stall-induced  oscillations  while  operating 
In  the  static  thrust  condition  showed  the  Initiation  of  a  pressure  disturbance  In  the 
region  of  the  leading  edge  as  the  blade  section  approaches  maximum  angle  of  attack,  and 
the  subsequent  motion  of  this  disturbance  In  the  chordwlse  direction,  at  considerably  less 
than  free-stream  velocity.  The  character  of  the  dlstu-bance  suggested  that  It  consisted 
of  free  vortlclty  Introduced  Into  the  blade  flow  flelo  from  the  neighborhood  of  the  blade 
leading  edge  during  the  dynamic  stall  process.  The  results  Indicated  that  the  dynamic 
stall  phenomenon  has  far  different  characteristics  than  those  associated  with  the  static 
stall  of  an  airfoil. 

The  negative  pressure  peak  generated  by  the  pressure  disturbance  moving  aft  from 
the  leading  edge  leads  to  a  nose-down  pitching  moment  component  In  phase  with  the  nose- 
down  motion  of  the  airfoil.  Since  this  nose-down  moment  Is  generated  once  per  pitching 
cycle,  It  is  seen  that  the  nonlinear  aerodynamic  moment  variation  due  to  pitching  motion 
at  high  mean  angles  of  attack  Is  such  as  to  sustain  the  motion.  This  self-excited,  self- 
limiting  motion  Is  termed  "stall  flutter". 

The  above  results  suggested  that  the  same  stall  mechanism  would  be  found  in  the 
general  case  of  large  transient  blade  angle-of-attack  changes.  Accordingly,  an  experi¬ 
mental  Investigation  was  commenced  at  MIT  to  study  large  linear  angle-of-attack  changes 
of  a  two-dimensional  wing  (Reference  (22)).  Comparison  of  the  dynamic  lift  variation 
with  the  corresponding  values  of  static  lift  at  the  same  angles  of  attack  Indicated  that 
the  maximum  dynamic  lift  achieved  was  substantially  higher  than  the  maximum  static  lift. 

A  very  large,  sustained,  transient  nose-down  moment  occurred  In  the  high  angle-of-attack 
region. 


The  origin  of  these  effects  was  found  In  the  corresponding  chordwlse  pressure 
variations.  Dynamic  stall  began  to  occur,  as  Indicated  by  the  drop  In  leading-edge  suc¬ 
tion,  at  an  angle  of  attac*  much  hlgner  than  that  for  static  stall.  A  negative  pressure 
disturbance  moving  aft  from  tne  leading  edge  simultaneously  Increased  the  suction  In  the 
mid-chord  region.  Subsequently,  the  pressure  disturbance  moved  further  aft  and  was  still 
of  considerable  magnitude.  The  delay  in  the  occurrence  of  stall  (as  evidenced  by  loss  of 
leading-edge  suction)  due  to  the  high  rate  of  cnange  of  angle  of  attack,  and  the  sustained 
upper  surface  suction  associated  with  tne  chordwlse  passage  of  the  vortlclty  shed  during 
tne  stall  process,  both  contributed  to  the  high  sustained  lift.  In  addition,  the  In¬ 
creasingly  aft  center  of  pressure  due  to  the  aft  motion  of  the  shed  vortlclty  generated 
extreme  nose-down  pitching  moment.  Finally,  the  pressure  distribution  Indicated  greatly 
increased  pressure  drag  on  the  airfoil.  This  drag  may  be  a  transient  analog  of  the 
"vortex  drag"  due  to  tne  leading  edge  vortex  of  a  slender  delta  wing  at  low  speed  and 
high  angle  of  attack. 

The  above  results  led  to  several  Important  conclusions  with  respect  to  3tall 
flutter  and  airload  prediction  of  high  speed  and/or  highly  loaded  helicopter  rotor  blades. 

1.  The  stall  of  an  airfoil  section  during  rapid  transient  high  angle-of- 
attack  changes  !s  delayed  well  above  the  static  stall  angle  and  re¬ 
sults  lr.  a  large  transient  negative  pressure  disturbance  leading  to 
large  transient  lift  and  nose-down  pitching  moment. 

2.  The  magnitude  of  the  pitching  moment  of  (1)  is  such  as  to  generate 
substantial  nose-down  pitching  displacements  of  the  blade.  These 
pitching  displacements  can  substantially  alter  the  angle-of-attack 
distribution  of  the  rotor  blade.  Transient  pitching  displacement 

of  the  blade  in  response  to  the  Initial  stall-induced  pitching  moment 
acting  on  the  blade  should  be  Included  In  stall  flutter  analyses. 

3.  The  dynamic  stall  phenomenon  of  a  helicopter  rotor  blade  can  be 
separated  Into  three  major  phases: 

a.  A  delay  In  the  loss  of  blade  leading-edge  suction  to  an  angle 
of  attack  above  the  static  stall  angle,  with  associated  airloads 
of  the  type  descrloed  by  classical  unsteady  airfoil  theory. 

b.  A  subsequent  loss  of  leading-edge  suction  accompanied  by  the 
formation  of  a  large  negative  pressure  disturbance  (due  to  the 
shedding  of  vortlclty  from  the  vicinity  of  the  blade  leading 
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edge)  which  moves  aft  over  the  upper  surface  of  the  blade. 

Associated  with  this  phase  are  high  transient  lift,  drag,  and 
nose-down  pitching  moment  associated  with  the  greatly  altered 
pressure  distribution  on  the  airfoil, 
c.  Complete  upper  surface  separation  of  the  classic  static  type, 
characterized  by  low  lift,  high  drag,  and  moderate  nose-down 
pitching  moment. 

b.  Method  of  Analysis 

Investigations  of  harmonically  oscillating  two-dimensional  wings  In  forced  motion, 
for  example.  References  (15),  (lb),  (17),  (18),  (19),  (20),  and  (21)  have  demonstrated  that 
under  stalled  conditions  the  average  damping  In  pitch  over  a  cycle  can  become  substantially 
negative  and  Is  strongly  dependent  on  the  wing  mean  angle  of  attack,  the  reduced  frequency 
of  the  harmonic  motion,  the  oscillation  amplitude,  and  the  airfoil  configuration  and  pitch- 
axis  location.  Reference  (16)  indicated  that  the  origin  of  the  negative  damping  was  aero¬ 
dynamic  moment  hysteresis,  and  that  for  certain  mean  angles  of  attack,  reduced  frequencies 
and  amplitudes  of  oscillation,  the  mean  damping  In  pitch  was  zero  over  a  cycle,  and  that, 
under  these  conditions,  a  self-excited  but  -self-limiting  one  degree  of  freedom  limit  cycle 
oscillation  of  prescribed  amplitude  could  occur.  These  results  Indicate  that  potential 
theory  unsteady  aerodynamic  predictions  of  two-dimensional  airfoil  pitch  damping  are  pro¬ 
gressively  less  representative  as  the  mean  angle  of  attack  approaches  the  static  stalling 
ang^e.  Typical  variation  of  the  pitch  damping  with  mean  angle  of  attack  Is  shown  In 
Figure  (11).  This  figure  represents  an  analytic  synthesis  of  data  presented  In  References 
(16),  (17),  and  (20).  In  the  referenced  experiments,  the  amplitudes  of  the  stable  limit 
cycle  oscillation  were  developed  as  functions  of  the  Initial  angle  of  attack  and  reduced 
frequency.  In  the  present  analytic  synthesis,  an  attempt  was  made  to  correct  for  differ¬ 
ences  In  rotation  point  and  static  stalling  angles  of  the  various  data.  This  was  followed 
by  a  calculation  of  equivalent  viscous  damping.  The  necessary  balance  of  energy  over  a 
cycle  of  oscillation  yielded  a  generalized  equivalent  viscous  pitch  damping  function. 

Thl3  equivalent  viscous  damping  function  by  the  nature  of  the  averaging  process  smooths 
out  the  higher  frequency  components  of  the  destabilizing  aerodynamic  pitching  moment.  It 
Is  Important  to  note  that  the  Instantaneous  values  of  negative  damping  vary  about  this 
mean  value  and  can  be  substantially  more  negative  near  the  stalling  angle. 

The  extent  of  the  stalled  regions  of  a  helicopter  rotor  and  the  possibilities  for 
unstable  pltchlng-torslonal  oscillations  are  Illustrated  in  F.  ;ure  (12)  which  Illustrates 
a  typical  angle-of-attacx  distribution  In  the  high  speed  cruise  condition.  This  distribu¬ 
tion  Is  based  on  a  detailed  calculation  of  the  rotor  velocity  field,  and  Indicates  an  ex¬ 
tensive  region  of  stall;  the  corresponding  net  Integrated  pitch  damping  Is  found  to  be 
negative  over  a  significant  range  of  azimuth  angles.  This  Indicates  the  origin  of  a 
transiently  unstable  pltchlng-torslonal  oscillation  which  will  occur  with  a  once  per  rotor 
revolution  repetition  rate.  Since  only  a  few  cycles  of  torsional  motion  are  possible  be¬ 
fore  the  blade  becomes  unstalled,  limit  cycle  motion  Is  usually  not  achieved.  However, 
substantial  Increases  In  blade  torsional  stress  and  pitch  link  loads  are  possible  either 
by  self-excltatlon  or  In  response  to  external  disturbances. 

Experimental  evidence  of  sucn  unstable  pltchlng-torslonal  oscillations  or  rotor 
blade  "stall  flutter"  Is  abundant.  The  term  "comfort  stall"  has  been  used  to  describe  an 
almost  asymptotic  rise  In  cyclic  pitch  link  loadings  and  related  helicopter  vibratory 
phenomena  which  occur  when  significant  zones  of  blade  stall  are  present.  Some  typical  ex¬ 
perimental  evidence  of  this  asymptotic  rise  in  cyclic  pitch  link  loadings  Is  presented  In 
Figure  (13)  In  the  form  of  a  set  of  wave  forms  showing  tne  tine  variation  of  the  torsional 
loading  In  a  typical  case. 

In  view  of  the  experimentally  derived  damping  function,  the  current  Insight  Into 
the  rotor  stalling  pattern,  and  this  experimental  evidence,  it  Is  evident  that  the  net 
pitch  damping  of  the  rotor  blade  can  and  often  does  become  periodically  negative  In  for¬ 
ward  flight.  Therefore,  whenever  combinations  of  thrust  ratio  and  rotor  advance  ratio 
result  In  significant  zones  of  blade  stall,  an  unstable  torsional  oscillation  may  result. 

The  stability  boundary  can  be  approximated  by  considering  the  net  aerodynamic 
pitch  damping  of  the  rotor  blade  control  system  fundamental  pltchlng-torslon  mode  of  os¬ 
cillation.  because  of  the  complex  angle-of-attacx  distributions  which  occur  In  forward 
flight,  this  net  damping  function  will  vary  widely  with  azimuth.  It  can  be  expected  to 
exhibit  a  heavily  damped  condition  In  the  region  of  blade  advance,  and  can  be  expected 
to  exhibit  a  negatively  damped  condition  In  the  region  of  blade  retreat.  If  large  values 
of  rotor  thrust  coefflclent-solldlty  ratio  and/or  rotor  advance  ratio  lead  to  significant 
zones  of  stalling.  Inasmuch  as  the  Instantaneous  negative  pitch  damping  can  exceed  the 
averages  shown  In  Figure  (11),  a  simple  approximation  of  the  stability  boundary  Is  ob¬ 
tained  by  the  condition  that  the  motion  will  be  transiently  unstable  If 
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where  the  reduced  frequency  Is  based  on  the  actual  velocity  at  the  reference  radius  of  the 
retreating  blade  rather  than  on  the  mean  velocity,  1 .  e . ; 
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In  other  words,  the  local  blade  element  average  damping  ratio  (depending  on  the 
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local  Initial  angle  of  attack  and  local  reduced  frequency)  Is  weighted  by  the  square  of 
the  local  fundamental  mode  amplitude  to  obtain  the  net  pitch  damping  at  each  azimuth 
angle.  If  this  weighted  average  damping  becomes  negative  at  any  azimuth  angle,  the 
pltchlng-torslonal  motion  can  be  expected  to  become  transiently  unstable.  If  the  range 
of  azimuth  angles  over  which  the  pitch  damping  Is  negative  Is  broad  enough  to  permit  one 
or  more  cycles  of  a  torsional  oscillation,  a  marked  increase  In  cyclic  control  loading 
can  be  expected. 

To  Illustrate  the  application  of  this  stability  criterion,  a  calculation  of  the 
net  damping  function  versus  azimuth  angle  was  carried  out  In  Reference  (20)  for  a  severe 
rotor  loading  condition  characterized  by  u  -  .17  and  CT/o  -  .111.  The  result  of  this 

calculation  Is  presented  In  dimensional  form  lr.  Figure  (1A).  It  Is  seen  that  the  net 
pitch  damping  Is  negative  for  the  region  bounded  by  the  azimuth  angles  225°  and  10°  (or 
370°).  Stall  flutter  would  be  expected  to  occur.  Figure  (13)  presents  torsional  strain 
and  pressure  data  for  this  flight  condition  which  Is  seen  to  support  the  theoretical  pre¬ 
diction  of  stall  flutter.  The  trace  of  absolute  pressure  transducer  at  the  80  percent 
radius  station  and  5  percent  chord  point  Is  Indicative  of  the  loss  of  leading  edge  suction 
and  accompanying  pitching  moment  variation.  It  Is  seen  that  the  blade  torsional  response 
to  this  Initially  nose-down  moment  exhibits  an  unstable  behavior  In  the  region  In  which 
the  net  damping  In  pitch  Is  negative.  In  this  case,  the  Instability  Is  short-lived,  but 
especially  pronounced,  as  evidenced  by  the  peak  In  the  pitch  link  load  trace  near  the 
330°  azimuth. 

A  more  precise  method  of  evaluating  stall  flutter  effects  has  recently  been  de¬ 
veloped  (References  (23),  (2A),  (25),  (26)).  The  equations  of  motion,  e.g.,  Subsection  2a, 
are  solved  numerically  for  the  actual  blade  motion.  Including  the  forcing  aerodynamic  terms 
due  to  blade  dynamic  stall. 

The  representation  of  blade  dynamic  stall  Is  based  on  the  results  of  References 
(22)  and  (27),  some  of  which  are  presented  In  Figure  (15).  The  applicability  of  these 
two-dimensional  results  to  the  rotating  blade  is  demonstrated  In  Reference  (28).  It  Is 
seen  that  under  conditions  of  rapid  transient .  as  opposed  to  oscillatory .  blade  angle-of- 
attack  change,  the  maximum  lift  and  torsional  moment  generated  are  proportional  to 

a^c/V,  where  (s)  denotes  the  value  at  the  Instant  of  dynamic  stall,  and  are  considerably 
larger  In  magnitude  than  those  measured  In  tests  of  airfoils  oscillating  through  the  stall. 
In  the  latter  case  It  Is  believed  that  some  degree  of  flow  separation  persists  throughout 
the  oscillatory  cycle,  leading  to  forces  and  moments  that  are  lower  than  those  found  In 
the  transient  tests. 

The  three  degrees  of  freedom  considered  In  the  analysis  are  rigid  blade  pitching 
about  the  feathering  axis,  rigid  blade  flapping  about  the  zero  offset  flapping  hinge,  and 
olade  first  mode  flatwise  elastic  bending.  Variable  Inflow  and  reverse  flow  are  Included 
In  the  computation  of  blade  forces  and  momenta.  The  calculation  of  Inflow,  aerodynamic 
loading,  and  blade  motion  Is  performed  Iteratively,  using  estimated  steady-state  Initial 
values . 

A  typical  result  Is  shown  In  Plgure  (16).  The  corresponding  rotor  angle-of- 
attack  distribution  Is  shown  In  Figure  (17).  Note  that  the  blade  passes  In  and  out  of 
stall  several  times  as  a  result  of  its  torsional  response  to  the  Initial  dynamic  stall. 

c.  Mear^  of  Alleviating  Stall  Flutter 

Several  approaches  to  the  design  problem  of  postponing  rotor  blade  stall  flutter 
are  evident.  The  first  Is  the  obvious  one  of  avoiding  significant  regions  of  blade  stall 
by  Increasing  rotor  solidity  beyond  that  dictated  by  ordinary  performance  considerations. 
This  approach  may  Incur  a  performance  penalty.  Excess  solidity  Increases  rotor  profile 
drag,  resulting  In  a  reduced  rotor  lift  to  effective  drag  ratio  at  the  design  point,  with 
an  attendant  loss  In  the  helicopter's  range  and  maximum  speed  capability.  A  second  ap¬ 
proach  Is  to  counteract  the  negative  aerodynamic  pitch  damping  with  positive  mechanical 
damping.  The  difficulty  stems  from  the  deformation  pattern  associated  with  the  funda¬ 
mental  torsion  mode.  Generally,  this  mode  Is  twisting  motion  In  addition  to  rigid  body 
pitching  against  the  effective  torsional  spring  of  the  helicopter  control  system.  In  this 
case,  a  dashpot  In  parallel  with  the  control  system  Is  relatively  ineffective,  and  In¬ 
ternal  torsional  damping  becomes  necessary;  this  Is  slight  In  conventional  metal  rotor 
blade  structures.  On  the  other  hand,  the  use  of  a  glass  fiber-resin  matrix  structure 
could  be  expected  to  dissipate  considerable  energy  In  torsion  due  to  its  visco-elastic 
damping  which  varies  with  the  rate  of  sheer  strain. 

The  most  appealing  approach,  wlti  the  least  performance  penalty  and  design  com¬ 
plication  is  to  reduce  the  extent  of  the  stalled  zones  by  departing  from  the  typical  con¬ 
temporary  blade  airfoil  section,  and  employing  sections  having  high  dynamic  stall  angles. 
Sections  having  an  appreciable  amount  of  leading  edge  camber  have  favorable  dynamic  stall 
characteristics  (References  (29)  and  (30)). 
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Fig.3  Effect  of  advance  ratio  /a  on  unsteady  aerodynamic  effects 
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Fig. 4  Flutter  and  divergence  boundaries  for  Happing  rotor  blade 
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Fig.6  Blade  response  to  classical  flutter  n  =  .68,  V  =  284  knots,  0O  =  4°  blade  C.G.  at  30%  chord 
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Fig.7  Blade  response  to  torsional  divergence,  p  =  1.47,  V  =  332  knots  0P  =  2°  blade  C.G.  at  257f  chord 
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Fig. 8  Geometry  of  undeformed  blade  and  deformed  blade 
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Fig.  I  2  Theoretical  angle  of  attack  distribution  with  non-uniform  downwash;  Vertol  model  CH-47A  front  rotor 
V  =  140  kts;  forward  C.G.  position,  gross  weight  =  27.500  lbs. 
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